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Abstract. The superconducting-insulator transition is simulated in disordered networks of Josephson junc- 
tions with thermally activated Arrhenius-like resistive shunt. By solving the conductance matrix of the 
network, the transition is reproduced in different experimental conditions by tuning thickness, charge den- 
sity and disorder degree. In particular, on increasing fluctuations of the parameters entering the Josephson 
coupling and the Coulomb energy of the junctions, the transition occurs for decreasing values of the crit- 
ical temperature Tc and increasing values of the activation temperature To. The results of the simulation 
compare well with recent experiments where the mesoscopic fluctuations of the phase have been suggested 
as the mechanism underlying the phenomenon of emergent granularity in otherwise homogeneous films. 
The proposed approach is compared with the results obtained on TiN films and nanopatterned arrays of 
weak-links, where the superconductor-insulator transition is directly stimulated. 



Although superconductivity has been predicted to per- 
sist even in the presence of impurities [1 , a supercon- 
ducting insulating transition (SIT) has been observed in 
strongly disordered metals |2|3|4|5|6|7|8j . A 'fermionic' phe- 
nomenon explaining the superconducting-insulating tran- 
sition is related to the suppression of the Cooper pair- 
ing by the enhancement of the Coulomb repulsion be- 
tween electrons with increasing disorder. An alternative 
'bosonic' phenomenon driving the transition is related to 
the phase fluctuations that could cause pair breaking and 
destroy the superconducting energy gap. Spatial inhomo- 
geneities of the order parameter might lead to complex 
effects due to the formation of islands with nonnegligible 
superconducting order parameter embedded in an insulat- 
ing matrix J9I1QI11I12I13I14J . 

Recent investigations carried out in homogeneous thin 
films by tuning disorder, electric or magnetic field pointed 
to a mechanism of self-induced granularity [15 16 17 18 19 . 
The emergent collective insulating state exhibits threshold 
voltage depinning behavior and thermally activated resis- 
tance described by an Arrhenius-like exponential law. In 
spite of the homogeneity, such systems have been thus 
envisioned to behave as networks of Josephson junctions, 
where the onset of the insulating phase requires the Coulomb 
energy Ec: 



En = 



2C 



(1) 



be larger than the coupling energy of the Josephson junc- 
tion: 

Different mechanisms leading to the superconducting-insulating 
transition very likely coexist in such complex electronic 
systems with collective insulating states emerging within 
the homogeneous structure. 

In this work, the "self-induced granularity" is modelled 
by disordered networks of Josephson junctions [20 21 22 2 3|24|25|26|' 
The superconductor-insulator transition is simulated by 
solving the conductance matrix of Josephson junctions. 
The networks are biased by a constant current generator 
/, as sketched in Fig. [T] The superconductive, intermedi- 
ate and resistive state of each Josephson junction is set ac- 
cording to each single current- volt age characteristic which 
is shown in Fig. [2] The approach was previously used to 
investigate the normal-superconductive transition and the 
current fluctuation power spectra in [28 . 

It is shown that an Arrhenius-like dissipation mech- 
anism in the shunted Josephson junction induces the in- 
sulating phase. Furthermore, the resistance variation ob- 
served in the array well compares with the experimental 
results of the superconducting-insulator transition in ho- 
mogeneous and nanopatterned TiN films [ 17fl8j . As al- 
ready mentioned, the resistance in the insulating state de- 
pends on temperature according to an Arrhenius-like ex- 
ponential as found in many experiments [15|16|17|18|19j . 
The resistance R of the Josephson junction array is pro- 
portional to the resistance Rj of each weakly linked Joseph- 
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Fig. 1. Scheme of a three-dimensional network of Joseph- 
son junctions. Circles are the superconducting nodes. The 
links between nodes behave as Josephson junctions. A^^^, 
Ny and A^;^ indicate the number of Josephson junctions 
respectively in the x, y and z directions. The network is 
biased by a constant current source with intensity /. 
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Fig. 2. Scheme of the current- voltage characteristic of 
a Josephson junction for the superconductor-normal (up- 
per branch) and insulator-normal (lower branch) states. 
The dashed line is obtained when the external drive de- 
creases as an effect of the hysteresis. Due to the disorder 
of the network, a wide variability of the I-V characteristics 
should be assumed according to the probability distribu- 
tion function of the relevant parameters of the Josephson 
junctions. In the present work, a Gaussian probability dis- 
tribution has been considered. 



son junction according to the relationship: 

N,N, + 1 



R 



Rj- 



N,, 



(3) 



where the subscript j identifies a single junction, while 
Nx, Ny and Nz indicate the total number of Josephson 
junctions respectively in the x, y and z directions as shown 
in Fig. IT] The resistance Rj of each Josephson junction is 
taken oithe form: 



Rj = Roj exp 



T 



(4) 



with Toj the activation temperature and Rqj the normal 
state resistance of each junction. The exponent 7 takes val- 



ues in the range to 1. Specifically, 7 = 1 corresponds to 
the Arrhenius law, while fractional values of 7 feature vari- 
able range hopping transport. Experiments have shown 
that 7 = 1 holds on the insulating side of the transition 
at low temperatures, while variable range hopping trans- 
port with exponent 7 = 1/2 or 7 = 1/4 are likely to occur 
over the normal side of the transition at higher temper- 
atures. The activation temperature Tqj is related to the 
parameters of the junction through: 



■Ecj/4. 



(5) 



This relationship means that a Cooper pair, in the su- 
perconducting islands separated by insulating barrier of 
each junction j, is firstly broken by supplying an energy 
Aj. Then the charge is redistributed in each tunnel ca- 
pacitance Cj by supplying an additional amount of en- 
ergy Ecj/4:. When the Coulomb energy is larger than the 
Josephson coupling, the value of the activation tempera- 
ture is thus mostly determined by the value of the junction 
capacitance. 

It is worth noting that the dimensionless tunneling 
conductance g of the Josephson junction is defined as: 



G 



(6) 



where G is the average conductance between adjacent 
grains and e'^/h is the quantum conductance. The aver- 
age conductance G of the array can be obtained as the 
reciprocal of average resistance R defined by Eq. Q. 

It has been observed both theoretically and experi- 
mentally that when the average tunneling conductance G 
is greater than the quantum conductance, i.e. when ^ ^ 1, 
which corresponds to strong coupling between grains, a su- 
perconducting state is achieved at low temperature. Con- 
versely, when ^ <C 1, which corresponds to weak coupling 
between grains, an insulating state emerges at low tem- 
perature [26]. This phenomenon can be explained by the 
existence of an additional dissipative channel due to single 
electrons tunneling between grains, resulting in an effec- 
tive Coulomb energy: 



E, 



c 



A 



(7) 



By comparing Eqs. ([2| and ([7|), one can note that Ej is al- 
ways larger than Ec for g ^ 1^ implying a superconduct- 
ing ground state. Conversely, for ^ <C l^ Ej is smaller than 
Ec, allowing the onset of an insulating state. Both the 
numerical and qualitative description clearly show, that 
a superconducting insulator transition occurs as Ej com- 
pares with ksT (disorder and temperature induced su- 
perconducting insulator transition) which ultimately cor- 
responds to compare Josephson and Coulomb coupling 
through the average tunneling conductance G of the array 
and the consequent coexistence of the fundamental mecha- 
nisms underlying the superconducting insulator transition 
described above. 

The state (superconductive, normal, insulating) of each 
Josephson junction in the array is set according to each 
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current voltage (I-V) characteristics. In the I-V curve rep- 
resented in Fig. l2| Icj and Vcj = IcjRj represent the 
critical current ana the critical voltage of a single Joseph- 
son junction j. For the superconductive branch (upper 
branch) of the I-V characteristics, the normal resistive 
state is achieved when the bias current Ij flowing through 
the junction j exceeds the critical current Icj. For the in- 
sulating branch (lower branch) of the I-V characteristics, 
the normal resistive state is achieved when the bias volt- 
age Vj across each junction j exceeds the critical voltage 

Vcj- 

The critical current is a function of temperature ac- 
cording to the following linearized equation: 



T 



(8) 



where Icoj is the lowest temperature critical current. An 
analogous relationship holds for the critical field Hcj. To 
simulate the network disorder, the critical current Icj^ the 
critical field Hcj and the activation temperature Toj are 
taken as random variables, which fiuctuate around their 
average value according to a Gaussian distribution with 
standard deviation a [28 . 

The average resistance R of the network depends on 
elementary shunt resistance Rj through the linear rela- 
tionship Eq. (Is]). Therefore, the variation AR^ which is a 
measure of the slope of the transition, depends on bias cur- 
rent /, on temperature T and magnetic field H^ through 
the variation of the elementary shunt resistance ARj of 
each weak-link: 



AR, 






H- 



dT 



(9) 



The terms on the right hand side of Eq. ([9| will be now 
related to the relevant parameters of the junctions Icj 
and Hcj used in the simulation. The first and second 
terms on the right hand side of Eq. ([9| can be written 
as {dRj/dQAIc and {dRj/dHc)AHc. The term dRj/dT 
can be written as derivative of a compound function: 



dRj _ dRj die 



dRj dHc 



(10) 



Moreover, from Eq. ([8|, dIc/dT c^ -Ico/Tc. A similar 
relation holds for the derivative of the critical field with 
respect to T. On account of the relation Ale = Ic—Ico oc a 
and AHc = He — Hco oc a, one obtains that the slope of 
the resistive transition AR is related to the disorder degree 
a. Since it is {dRj/dQ < and {dRj /dHc) < 0, Eq. ^ 
shows that if either Ale or AHc increases, ARj decreases. 
The quantity ARj turns out to be proportional to the 
slope of the resistive transition curve, which is smoother 
when either Ale or AHc increases, i.e. when the disor- 
der increases. Therefore, the disorder enters the simula- 
tion directly by varying the parameter a. The disorder de- 
gree is strongly dependent upon the temperature T along 
the transition as the Josephson phase fiuctuations depend 
on T. For the sake of simplicity and because the work is 



mostly aimed at reproducing the transition from the su- 
perconducting to the insulating phase at a given low tem- 
perature, the simulations have been performed at constant 
value of a along the superconducting normal transition 
curve. 

Additionally, the SIT transition is shown to occur ei- 
ther when thickness decreases or when current increases, 
while the parameter a is kept constant. The aim of these 
simulations is to show that external drives, such as bias 
current or thickness variation, directly act on the network 
as "disorder enhancer" and trigger the onset of the insu- 
lating phase. The temperature dependence of the disorder 
{a) would further accelerate the occurrence of the insulat- 
ing phase in such cases as well. 

The simulations of the superconductor-insulator tran- 
sition in the Josephson junction arrays are performed as 
a function of the temperature T in the range to IK by 
using different tuning parameters, as electric field, thick- 
ness and disorder degree. The resistive shunt is taken of 
the standard Arrhenius law form with 7 = 1. 

Figure [3] shows the superconductor-insulator transition 
when the temperature increases and the bias current / 
is kept constant in two-dimensional disordered networks. 
Different curves correspond to different disorder, i.e. to 
different values of a. By increasing temperature, the cur- 
rent Icj of each junction decreases according to Eq. Q. As 
Icj becomes smaller than the bias current fiowing through 
each junction, the junction becomes resistive. In the ideal 
case of a network with no disorder, it is a = and the tran- 
sition curve is vertical. All the Josephson junctions have 
the same critical current, become resistive at the same 
temperature and the transition occurs simultaneously all 
through the network. For the curves shown in Fig. [3J the 
average value of the normal resistance Rq is taken equal 
to 20kl] and the average value of the critical temperature 
Tc = 0.8 for a = 0. Moreover, when <j = the activated 
conduction processes are negligible, so that the average 
value of the activation temperature can be taken equal to 
zero. To — 0. The average value of the critical current Ic^o 
is 20/iA. 

As disorder increases {a increases accordingly) , the re- 
sistance changes more smoothly due to the wider spread of 
the Icj values. Upon further disorder increase, the transi- 
tion becomes smoother and smoother and for even larger 
disorder the network becomes an insulator. As observed in 
many experiments, in the case of the standard Arrhenius 
law, the average value of the activation temperature To is 
of the same order of magnitude of Tc. Both Tc and To are 
related to the disorder in the film, however Tc decreases 
while the average value of To increases as disorder in- 
creases. One can assume To oc a implying that To increases 
when disorder increases. By increasing a the value of the 
critical temperature Tc decreases accordingly. The average 
values of the activation temperature To corresponding to 
each curve range from 0, for a = 0, to 8, for a = 1, with 
step 1. The values of the parameters used in these simula- 
tions are in the range of experimental values reported in 
il7|18j. 
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Fig. 3. Superconducting-insulator transition of a two- 
dimensional network of Josephson junctions as a function 
of temperature. For all the curves, it is N^ = 30, Ny = 50, 
A^;^ = 1. The bias current is kept constant / = 9.5/iA. 
Different curves correspond to different degrees of disor- 
der with the standard deviation of the critical currents a 
varying from to 1 with step 0.1. 
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Fig. 4. Superconducting insulator transition of a two- 
dimensional network of Josephson junctions as a function 
of temperature. For all the curves, it is N^ = 30, Ny = 50, 
Nz = I. The degree of disorder is kept constant by taking 
the standard deviation of the critical currents a equal to 
0.3. Different curves refer to different values of the bias 
current / ranging from 11/iA to 22/iA with steps 1/iA. 



Figure |4] shows the superconducting insulator transi- 
tion as a function of the temperature, for different val- 
ues of the bias current /. The degree of disorder is kept 
constant by taking constant the value of a (in particular 
(J = 0.3). By increasing the bias current, the resistive tran- 
sition becomes smoother and the superconductor insulator 
transition occurs at larger values of the bias current. 

Figures |5J [6] and [7| show the superconductor insulator 
transition as a function of the temperature for different 
values of the film thickness, which is varied by changing 
the number of junctions Nz in the z-direction. The bias 
current is constant (/ = 50/iA). In each figure, the degree 
of disorder is constant, respectively by a = 0.3, a = 0.2 
and a = 0.1. One can notice that the simulations are ro- 
bust with respect to the variation of Nz at different disor- 
der degrees and the insulating phase appears for thinner 
films (i.e. for thinner networks corresponding to smaller 
values of TV;^). 

In this paper, the superconducting-insulator transition 
has been simulated in arrays of thermally activated resis- 
tively shunted weak-links at varying disorder, bias current 
and film thickness. Accurate predictions of the relevant 
properties of the insulating state have been obtained. In 
particular, the agreement between simulation and exper- 
iment is consistent with a complex process according to 
which the homogeneous system self- organizes into a gran- 
ular structure of weak-links |15|16|17|18|19j . The onset of 
the insulating phase, whose relevant property is the ex- 
ponential dependence of the resistance R on temperature 
T, given by Eq. Q, is recovered, consistently with the 
existence of a nonvanishing gap in the insulating phase. 
Further applications and developments of this work can 
be envisioned to shed light on the elementary processes 
underlying the superconducting-insulator transition. The 
flexibility and ease of implementation of the proposed ap- 




T(K) 

Fig. 5. Superconducting-insulator transition of a three- 
dimensional network as a function of temperature. For the 
curves shown in this figure, A^^^ = 20, Ny = 25 whereas 
Nz is varied from 1 to 10 in steps of 1 to simulate the 
film thickness variation. The bias current is kept constant 
to the value / = 50/iA. The degree of disorder is kept 
constant with the standard deviation a = 0.3. 



proach would be interesting for further investigations of 
the emergence of such complex electronic structures in 
several frameworks [29.30,31 32j. 
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